Spectral properties of quasi-one-dimensional conductors with a 
finite transverse band dispersion 

Z Bonacic Losic 1 , A Bjelis 2 and P Zupanovic 1 
1 Department of Physics, Faculty of Natural Sciences, 
Mathematics and Kinesiology, University of Split, 
Teslina 12, 21000 Split, Croatia^ and 
2 Department of Physics, Faculty of Science, 
University of Zagreb, POB 162, 10001 Zagreb, Croatia^ 

Abstract 

We determine the one-particle spectral function and the corresponding derived quantities for 
the conducting chain lattice with the finite inter-chain hopping t± and the three-dimensional long- 
range Coulomb electron-electron interaction. The standard GoWo approximation is used. It is 
shown that, due to the optical character of the anisotropic plasmon dispersion caused by the finite 
t±, the low energy quasi-particle 5-peak appears in the spectral function in addition to the hump 
present at the energies of the order of plasmon energy. The particular attention is devoted to the 
continuous cross-over from the non- Fermi liquid to the Fermi liquid regime by increasing t±. It is 
shown that the spectral weight of the hump transfers to the quasi-particle as the optical gap in the 
plasmon dispersion increases together with t±, with the quasi-particle residuum Z behaving like 
— (lnt^) -1 in the limit t± — > 0. Our approach is appropriate for the wide range of energy scales 
given by the plasmon energy and the width of the conduction band, and is complementary to the 
Luttinger liquid techniques that are limited to the low energy regime close to the Fermi surface. 
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I. INTRODUCTION 



Recent ARPES measurements of photoemission spectra show that a series of quasi-one- 
dimensional conductors, in particular the acceptor-donor chain compound TTF-TCNQ 
and Bechgaard salts (TMTSF) 2 X with X = PF 6 , C10 4 , Re0 4 ,... [3|, U, [5j, have unusual prop- 
erties, clearly distinguishable from the spectra of standard three-dimensional conductors. 
Quasi-particle peaks in these compounds are absent, and the spectra are instead dominated 
by a wide feature spread across energy scales of the order of plasmon energies. Such data 
are in qualitative accordance with the conclusions of our recent calculation [6| for the spec- 
tral function of the one-dimensional electron band with the three dimensional long range 
Coulomb electron-electron interaction, obtained within the so-called GoWo approximation 
The physical origin of such behavior is the one-dimensionality of the electron band that 
causes an anisotropic acoustic plasmon dispersion. Since such dispersion spreads through 
the whole range of energies, from zero up to the plasmon energy Q p i, it introduces the wide 
feature into the spectral function at these energies, leaving thus no space for the creation of 
quasi-particle 5-peaks. 

The spectral density N(u) and other quantities related to the electron spectral properties 
have been also calculated exactly within the Luttinger liquid approach, using mostly the 
bosonization method 8, l|. Such analyzes are however limited to the narrow range of low 
energies, u <C Ep, Q p i, where Ep is the Fermi energy of the order of bandwidth. It was 
shown that, together with the absence of quasi-particle peaks, the spectral function shows 
power law behavior with the anomalous dimension a, defined by N(u) ~ \uj\ a and 
being interaction dependent. The comparison with measurements at low frequencies suggests 
values of anomalous dimension in the range a > 1 . This corresponds to the regime of strong 
three-dimensional long-range Coulomb interactions 1(3, 11, 12, Q], which additionally 
suggests that the corresponding plasmon energy scale is not small, being at least of the order 
of band-width or larger. The GqWq approximation is the only known approach which, as 
was already pointed out, enables the calculation of spectral properties in such wide ranges. 
However, it does not lead to the correct power law exponent in the limit ui — > 0. As such, it 
is complementary to the Luttinger liquid approach jg, which is concentrated and limited 
to the low energy region. 

The combination of two above approaches thus covers the whole energy range relevant 
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for the analysis of the photoemission properties of quasi-one-dimensional metals. As was 
already stated, the main emerging conclusion for the electron liquid with a strictly one- 
dimensional band dispersion is that, although three-dimensionally coupled through long- 
range Coulomb interaction, it does not show the essential property of Fermi liquids, namely 
the presence of quasi-particle excitations in the one-particle spectral properties. However in 
order to understand better the spectral properties of real quasi-one-dimensional conductors 
one has to take into account deviations from the one-dimensional band dispersion which come 
from finite inter-chain electron tunnellings. The corresponding question of both, theoretical 
and experimental interests is: how one reestablishes the Fermi liquid character of spectral 
properties by introducing and gradually increasing the transverse bandwidth t±, approaching 
thus the regime of standard isotropic three-dimensional conducting band? 

In this work we address this question by extending our earlier GqWq approach to the 
rectangular lattice of parallel chains with a finite transverse tunnelling integral After 



taking into account 



band dispersion [15 



;he corresponding finite transverse curvature in the three-dimensional 



161 ] . the screened Coulomb interaction Wq calculated within the random 



phase approximation (RPA) shows a finite optical plasmon gap proportional to t± in the 
long- wavelength limit. The plasmon dispersion thus has a three-dimensional, albeit strongly 
anisotropic, character for any finite value of t±. A more detailed insight into the electron 
self-energy within the GqWq approach shows that this property of plasmon dispersion has 
the dominant effect on the dressed electron propagator through the screened Coulomb in- 
teraction Wo, while the influence of finite t± through a bare electron propagator Go can be 
neglected. This enables an analytical derivation of the dressed electron Green's function and 
other quantities that follow from it. 

The obtained result reveals the appearance of low energy quasi-particle peaks, in addition 
to the smeared structure at higher energies which is a characteristic of the strictly one- 
dimensional (t± = 0) limit jo]. Note that the early GqWq approach to the isotropic three- 
dimensional ,,jellium" IT . Tisl. [3] led to the analogous result for the spectral function, 
showing quasi-particle peaks in the energy range fi — Q p i < uj < jj, + Q p i where fl p i is the 
minimum of the optical long- wavelength plasmon dispersion, and an additional structure 
due to the plasmon mode, with the finite spectral weight below and above these energies. 

The spectral properties for the generalized Luttinger liquid with a weak electron tun- 
nelling between metallic chains and with the three-dimensional electron-electron Coulomb 
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interaction were ana 
tion technique 
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yzed by using the appropriately developed higher-dimensional bosoniza- 



14j in which the Fermi surface is approximated by a finite number of 



flat patches. This technique inherits in itself two approximations, namely the momentum 
transfer between different patches is ignored and the local band dispersion is linearized. On 
the other hand, it handles the case of t± 7^ without having to rely on an expansion in 
powers of t \ t used in earlier studies of the model of parallel chains with a finite inter-chain 

241 ] . Using the 4-patch approximation for the Fermi surface Kopietz 




obtained in the strong coupling limit the spectral function with the low 
energy quasi-particle having the weight proportional to 6 7cb , 6 = \t±\/Ep. Here 7 c b is the 
anomalous dimension of corresponding Luttinger liquid for t± = 0, and Ep is Fermi energy. 
Furthermore, it is shown that there exists a large intermediate regime of wave vectors and 
frequencies where the Green's function satisfies the same anomalous scaling behavior as for 
t± = 0. This is to be contrasted with the result of the perturbation treatment of t± [2^] 
in which the quasi-particle peak appears only when the one-dimensional Green's function 
diverges, i. e. for the anomalous dimension less than unity. 

Again, like in the case t± = 0, the higher dimensional bosonization and our GqWq ap- 
proach are complementary, since the former is limited to the scaling behavior of the Green's 
function in the low energy range and the latter enables the reliable calculation of the wide 
maximum at the range of plasmon energy in the spectral function. It is important to note 
that the essential ingredient in both approaches is that the finite t± enters into calculations 
through the long wavelength optical gap in the plasmon dispersion, and not through the 
corrugation of the band dispersion at the Fermi energy as in the perturbation approach 
of Ref. 20]. On the other side, while both Wen's expansion in terms of t± 20j and the 



higher-dimensional bosonization treatment cover low energy scaling, only the present GqWq 
approach describes appropriately the cross over from the one- dimensional non-Fermi liquid 
regime to the three-dimensional Fermi liquid one in the whole range of energies. 

In Section II we calculate the electron Green's function within the GnWn method devel- 

n 

oped in our previous work [6|. Section III is devoted to the spectral function. The density 
and the momentum distribution function are discussed in Section IV. Section V contains 
concluding remarks. 
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II. GREEN'S FUNCTION 



A. Dielectric function and excitations 

We begin by considering the effect of finite transverse bandwidth on the plasmon disper- 
sion. The electron band dispersion is modeled by 

E(k) = — 2t (cos k\\b — cos k F b) — 2t^(cos k x a + cos k z c), (f ) 

where b and a, c are longitudinal and two transverse lattice constants respectively, while to 
and t± are corresponding transfer integrals. The RPA polarization diagram now reads 

Ufa co) = 4 V f f n(k)[jg(k + q)-2S(k)] 

iq ' ) N a N b N c k k^-* & =-tl i w + ^signa;) 2 — [E(\s. + q) — E(\s)) 2 ' [) 

x a II t) 2 c 

where 

n(k) = \ (3) 
[ 0, E(k) > E F 

is the occupation function. In the long wave-length limit q — > 0, where u; ^> J5(k+q) — E(\s), 



the polarization diagram reduces to 



with 



"> ^ w.(. 2 + w J, J, jL n(k)(q ' Vk)2B(k) (4) 



, „ , ,a 2 B(k) ,a 2 B(k) ,3 2 £(k) 



Since by assumption t± <C £o> the Fermi surface is only slightly corrugated, i. e. 
S(k x , k z )/k F <C 1, where 5(k x ,k z ) is the deviation of the component of the Fermi wave 
vector in the chain direction from k F , the latter being its value at t± = 0. The expansion of 



151 ] leads to the equation for 



the band dispersion (TjQ) in terms of 5 up to the second order 
the Fermi surface 

E(k x , k F + 5, k z ) = v F 5 + E F 5 2 /2 - 2t_L(cos k x a + cos k z c) = E F (6) 

where v F = 2to&sin k F b is Fermi velocity, E F = d 2 E(k)/dk 2 at k\\ = k F , and E F is the shift 
of the Fermi energy with respect to its value for t± = 0. Our aim is to find out how 5(k x , k z ) 
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depends on t±, and to determine the corresponding value of Ep. To this end we note that 
by switching to finite t± the band filling does not change, so that 

*7r/a rn/c 
r/a J—tt/c 

Then, since 5 ~ t± to the lowest order, the integration of Eq. [6] in terms of k x and k z gives 
Ep ~ t\. The explicit expansions follow after expressing 5(k x ,k z ) from Eq. [6j 



/n/a i"k/c 
/ 5(k X} k z )dk x dk z = 0. 
-n/a J—tt/c 



(7) 



v F 



El 



Ep — 2t^(cos k x a + cos k z c) 



2t± 

Vp 



(cos k x a + cos k z c) 



Ep 

Vp 



E 



— 2— |^(cos k x a + cos k z c)' 



(8) 



Inserting this expression into the condition ([7j) one gets Ep = 2E F t 2 ± /vp, and finally 
5(k x , k z 



2t i , . E p o E 'p o / i \ r 

(cosA^a + cosfc z c) + 2—^-t ± — 2— g-tj_(cos k x a + cosk z cy 



Vp 



(9) 



The expansion (Q enables the analytical derivation of the dielectric function e m (q, u) 



1 — y(q)Il(q, lu), where V(q) 



voq 2 



is the bare Coulomb interaction. After replacing 



-5*ii=-f ^ 



r \ 3 /.i 

lv \ /a 



dk 7 



k F +S 



dkn 



-(k F +S) 



in Eq. HI and taking into account that 



k F +8 



(k F +S) 



dk\ 



2 E(k) dE(k) 



dkl 



dh 



k F +S 



2(v F + E F 8 + E' F , 5 2 /2) 



(10) 



(11) 



and 



we get 



d 2 E(k) 2 , d 2 E(k) 2 , 

— — — — = 2t±a cos k x a, — — — — = 2t±c cos k z c, 
dkl dk 2 



^ 2 (q) 



e m (q,uj) = l-- . . 

(lo + irjsignuy 



with the plasmon dispersion given by 



^ 2 (q) = — ^ — • 



(12) 



(13) 



(14) 



Here longitudinal and transverse plasmon frequencies are given by ft 2 t = 8e a ^ F (l + 2-^f L tj_) 



and c<j„„ = — , lo* = — 

P a cvp ' V c avp 



respectively. Thus, the finiteness of the transverse bandwidth 



retains the optical character of the plasmon dispersion in all directions of the long wavelength 



6 




FIG. 1: Plasmon dispersion a>(q) (see Eq. [T5|) . 

range q — > 0, with the anisotropy scaled by the ratio ^ as shown in Fig. [1] As long as t± -C to 
we can skip the correction proportional to t\ in f^. Also, for simplicity we put a = c and 
get the simplified expression for the long wavelength plasmon dispersion, 

2 Kl4 + <4iqj 

u (q) = 5 ' t 15 J 



with uj^ = 16 ^ fj - and q\ = q 2 + q 2 z . Note that in the regime of strong Coulomb interaction, 
Q P i 3> [6;], we also have 

>1. (16) 



w P « _ fac Vl p i 



t± V 2b 2 t sm(kpb) 



B. Green's function 



In the calculation of the reciprocal Green's function G _1 (k, u), we follow the GqWq ap- 
proximation [6j. The extension of this procedure obtained by the inclusion of the full q- 
dependence in the band dispersion (JTJ) leads to the generalization of the equation (20) in 



Ref. 



G-\k, u)=u- E(k) + i V [l - 2n(k)] - E ex {k) - ^ £ V{q)u{q) x 



X 



1 — n(k + q) 



n(k + q) 



/i - oj(q) - E(k + q) + irj oj — ji + u;(q) - E(k + q) 



(17) 



Here 



E ex (k) 



N 



^l/(q)n(k + q) 



(18) 



is the exchange energy per elementary cell for the one particle state with the wave vector 
k. Further simplification follows after noticing that, as far as we are in the regime of strong 
Coulomb interaction, Q p i ^> to (see Ref. 6j and Eq. [T6i) . two second terms in the dispersion 
E(k + q) « -Eo(^n) + Vpq\\ + E±(k a + q a , k c + q c ) appearing in the denominators of Eq. [TT] 
can be neglected with respect to that of the plasmon dispersion cu(q). As it will be seen 
later, this approximation introduces small losses in the spectral density at low frequencies, 
but does not affect its main qualitative features. After a few nonessential simplifications 
which do not affect the physical content, like taking the flat Fermi surface at |faj| = kp for 
the occupation function ([3]) and using cylindrical coordinates in the integration across the I 
Brillouin zone {(J, one gets the analytical expression for G _1 (k, uS). Its real part reads 



iteG _1 (k,w) 



<" " E(k) + -{ In 



bQ 



J_ 



71 



+ 1 



2bQ ± vr 
H arctan 



7T 



bQ_ 



+ 



e 2 j (u - /i - E (k\\))u p i 2tt 
2ir\(u — n — E {k\\)) 2 - lo 2 x b 

- fl - Ep^)) 2 

(u — fi — E Q (k\\)) 2 - uo. 



In 



+ In 



P i I 



71 



,0J 



LUpl + Vlpl 

fij — R(k\\,u> — n) 



1 + 



bQ. 



Ti 



+ 



\ 



to 



P i 



Ti 



2tt 

+ — In 

b 



/i-So(fe|i)Wl + 



— a> 



n 2 



/i - -Eo(faj) - ^ 



pi 



a; 



//-Eo(fc||))n 



pi 



u pl ((u - v - E (k l{ )) 2 - n 2 ) J-^ 



bQ ± 



dy 



+Q_ 



{u- ii- E {k\\)f{uo 2 pl 



u p i{{u -ii- £ (faj)) 2 - Q 2 pl ) 2 J-^- y^qn 



7T 



dy 



n l'„2 7 1 u , (^-M-go(fcii)) 2 -^ 



(19) 



with functions i? and F given by the expressions 



R{k h u) 



Ri{kp — \k\\\,oj) + Ri(kF + |fai 



+ 



Ri(ki 



Ti 



8 (l 

2tt 



I fall — &7 



(20) 



fen |, cj) + 2i?i -,w — Rx I — fci? — 



7T 



9(^+1^1- -J, 



8 



with 



Ri(x, uj) 



and 



F(x,u) 



2Qj 



Q. 




x 2 + Q 



{u>-Eo(k\\))2-wf 



f- x arctan • 



(q>-g (fc||)) 2 -a>g, 

^ ; -(^-i?o(fc||)) 2 



x In 



x—Q A 



2 (^-i?o(fc||)) 2 -^ 
^( W -E (fc||)) 2 -^. 



+ x ^ 0, 



(21) 



'"- fi o( fc ll» 2 -"pi 

( W -B (fc,|))2_n2 



/or |a> - E (k\\)\<u p i, fl p i<\u - E (k\\)\, 





Bo( fc ||)) 2 " 


- 2 ) 




-(w-B Q (fc 


|)) 2 


/ (w- 


s C=||)) 2 - 


- 2 , 
pi 




-(w-E (fe 





/or L) p i<\u - E (k\\)\<Q p i. 



The exchange energy in the expression (1191) is given by 



(22) 



Eex(k\\) = ~ 



2tt 



H(k F -\h\\) + H(k F + \k\\\) 



e 



+ 



+ 



7T 



27T 



with 



#(fc F - |Jfeii|) + 2if - - H — - Jfe 



= xln(Qj_ + x 2 ) + 2Q± arctan 



7T 



e( a* + |*u | - - 



x\nx . 



(23) 



(24) 



The further simplification follows after realizing that in the regime of strong Coulomb interac- 
tion the self energy contribution is dominant in comparison to the transverse dispersion term 
2tj_(cosk x a + cosfc^c) in E(k). Consequently we can skip the dependence of -ReG _1 (k, uj) 
on k x and k z as irrelevant for further considerations. Namely, after taking into account 
that Q± = 2y / 7r/y / ac <C tt/6, the leading contribution to the third term on the right hand 
side in Eq. [19] reduces to ~ 0.16 ^'^ p ' <C t±. This justifies the above approximation, 

after which we can proceed to a great extend along the lines of Ref . 6] . In particular, the 
chemical potential fi in Eq. [19] is now, after taking into account the self-consistent condition 
ReG~ 1 (k F , /i) = 0, given by 



- — { In 

2b 



7T 



+ 1 



2bQ ± 7T 

H arctan ■ 



7T 



bQj 



(25) 



The imaginary part of the reciprocal Green's function is given by 



ImG' 1 



(u-fi- E (k\\))' 



2q c Q(u -pi- E (k\\))- e(-w + A* + ^o(*||)) 



2 (w-Ai-^W-wJ 
+G(u; — \i — E (k\\)) x 2q c Q(k F — |A:||| — g c ) + 2k F Q(q c — \k\\\ — fc_p) 
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(a) 



(b) 



(c) 



FIG. 2: Frequency dependence of ReG (k\i,uj)/to (full lines) and —ImG l {k\\, w)/to (dashed lines) 
for fcp = 7r/26 and kn = (a), kn = k F (b), and kn = 2k F (c). 



2tt 

+(k F - \k\\\ + q c )Q(\k\\\ +q c - k F )Q(k F - \ \k\\\ - q c l)©(-^- - k F - \k\\\ - q 

+{2k F + 2q c - —)Q(k F - | \k\\\ - q c |)6( — r + A;f + N + g c 
o o 



(26) 



for tUpi < \u — fi — E (k\\)\ < Q p i, and ImG 1 (k\\,uj) = elsewhere. The wave number q c in 
Eq. [26] is defined by 



q c = mini Qj_ 



\W pl -(u;-»-E (k ll )rb 



(27) 



iteC? (fty, to) and JmG (fen, cj) are shown in Fig. [2]for three representative values of kn, 
namely for kn equal to 0, k F , and 2k F . Let us at first look more closely into ImG^^kn, uj). 
The vanishing of ImG -1 (k \ \, uj) in the ranges \uj—fi — Eo(k\\)\ < u p i and \u — fj l —Eo(k\\)\ > Vt p i 
can be traced already from the expression (ITT)) after approximating i?(k + q) ~ E (k\\). 
Namely, in this ranges there are no poles of the reciprocal Green's function contributing to 
ImG- l {k h uj). 

ImG~ l (fey, uj) vanishes also in the range [i+uj p i+ -Eo(^||) < uj < fj,+uj(kn—k F ,Qx)+E (kn) 
for fcy < k F , as well as in the range \i — uj(k\\ — k F , Q±) + E (k\\) < uj < \i — uj p i + E (k\\) 
for fen > k F . This vanishing can be also traced from the expression (ITT)) . Namely, due to 
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the presence of the occupation function n(k + q) in the q - summation the non-vanishing 
contributions from dense discrete poles at uj = fi — uj(q) + Eo(k\\) + irj contribute only in 
the range /i — Q p i + E (k\\) < uj < ji — uj{k\\ — k F ,Qj_) + E (k\\), while the non-vanishing 
contributions from poles at uj — fi + uj{q) + E Q (k\\) — it] contribute only in the range ft + 
u(k\\ - k F , Q ± ) + E (k\\) < uj < fi + Q p i + E (k\\). 

In the range uj p \ < \uj — fi — Eo(k\\)\ < Q p i ImG~ 1 (k\\,uj) is covered by the expression 
( |26|) . It has a step singularity of the width Q - 2 %' at uj — fi ± Vt pl + E Q (k\\) and diverges 

pi pi 

at the energies uj = ji — uj p i + E Q (k\\) for fey < fc^ and uj = fi + uj pi + £?o(^||) f° r ^|| > 
&f- At energies = fi =F tj(ir/b, Q±) + E (k\\) ImG^ 1 (k\\,uj) has respective anomalous 
minimum and maximum, with jumps in the first derivatives. These extrema originate from 
the confinement of the q - summation in the expression ffTTl) to the first Brillouin zone. The 
integration in terms of q± from to Q± results in the limitation on the q\\ - integration 
to the range \q\\\ < Q L\j %- ^^°^ Fo ^y2 as far as this limit is within the I Brillouin zone. 
However, for values of uj in the ranges fi — Q p i + E (k\\) < uj < fi — u(7r/b, Q±) + -Eo(^||) and 
fi + uj(Tc/b,Q±) + E (k\\) <uj < fi + Vt p i + E (k\\) we have f < Q±yJ ^ ^f^E {k ^2 , so that 
the qu - integration is limited to the I Brillouin zone, i. e. by the uj - independent boundary 

q c = ?. The resulting values of ImG~ 1 {k\\, uj) at the anomalous minimum and maximum 

are Tf 2 L K.2-M-go(fc||)) 2 

Let us now consider ReG~ l (k\\, uj). As is seen from Fig. [21 it diverges towards ±00 at 
the respective energies uj — fi =p Q p i + E (k\\) at which ImG~ 1 (k\\, uj) has step singularities. 
These singularities are shifted towards larger values of uj as k\\ increases. The zeroes of 
i?eG~ 1 (A;||, uj) at uj < fi — Q p i + E {k\\) and uj > ft + Q p i + E (k\\) are also shifted to the 
right as k» increases, former approaching the singularity at uj = fi — Q p i + E (k\\) and latter 
increasing the distance from the singularity at uj = fi + Q p i + E (k\i). ReG~ l (k\\,uj) has also 
essential singularities at uj = fi — uj p \ + E (k\\) (for k\\ < kp) and uj = fi + uj p i + E Q {k\\) (for 
k\\ > kp), i. e. at energies at which ImG~ 1 {k\\, uj) diverges. 

The zero of ReG^^k^, uj) in the range ft — uj p i + E (kn) < uj < ft + uj p i + E (k\\) in which 
ImG~ 1 {k\\ 1 uj) vanishes is the low energy pole of the electron propagator G(h\,uj). It is of the 
form y{k\\) = E(k\\) — zT(A;||), where r(fc||) is infinitesimally small in the present approach. 
Accordingly, our Green's function has in this range the standard resonant form 
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FIG. 3: Zp obtained numerically (full curve) and from the expression (|29j) (dashed curve). 

where Z(k\\) —\ dReG~ l {k\\,y{k\\)) / du) \~ l is the residuum of the Green function at the pole 
y(k\\). We emphasize that the low energy pole appears due to the optical gap u p i in the long 
wavelength plasmon dispersion introduced by the finite interchain transfer integral t± in the 
electron dispersion. This is illustrated by the analytical expression for the residuum Z(k\\) 
at fen — kp in the limit uj p i -C fl p i, 

z r= i/h + £9±U^)] =1 /\ 1 + i.Q± bl (r$e&)] . (29) 

IV 7T \l p i V uj p i J \ IV 7T n p i \ et± J . 

The dependence of Zp on t± obtained numerically, as well as with the use the expression 
( I29p . is shown in Fig. [31 The Green's function has the standard resonant form ( |28l) also in the 
frequency range \u — (/, — Eo(k\\)\ > Q p i in which i?eG _1 (/c||, a;) has zeroes and JmG _1 (/c||, uj) 
vanishes. 

On the other hand the structure of the Green's function in the region uo p i < \u — fi — 
Eo(k\\)\ < Q p i in which ImG~ l (k\\, uj) 7^ is influenced by the plasmon dispersion contribu- 
tion to the expression (|T7|) . 
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III. SPECTRAL FUNCTION 



The single-particle spectral function is defined by 

A(k h uo) = - | ImG(k\\,uj) | . (30) 

It can be directly expressed in terms of i?eG _1 (fc||, uj) and imG _1 (/c||, uj), 

, . 1 | ImG^ 1 (k\\, uj) | 

( U) = n [ReG-\k h u)f + \ImG-\k h uj)Y' { } 

unless in the case of ReG~ 1 {k\\ 1 uj) having a zero y{k\\) in the frequency range in which 
ImG~ 1 (ku, uj) = 0, when it is represented by the quasi-particle 5-peak 

A{k h u)=Z(k\{)6(u-y(k\\)). (32) 



The spectral function A(k\\,uj), obtained after inserting expressions (jT9l) and (126|) into 
Eqs. fl3Tl) and (132]) . is shown in Fig. H]for two values of the transverse plasmon frequency, 
uj p i = 0.26eV and 0.63eV. Generally it is characterized by the coexistence of wide humps 
and quasi-particle 5-peaks. Humps originate from the plasmon dispersion in the range 
uj p i < \u — jj, — E (k\\)\ < Qpi. Their positions vary slowly with the wave number kn. As for 
the #-peaks, they are situated in the energy ranges /x + E (k\\) — u p i < uj < /i + E (k\\) + u p i 
and \u — \i — E (k\\)\ > Q p i. It is to be noted that 5-peaks are present for any finite t±. 
However, the decrease of t± leads to the decrease of the weight of the quasi-particle 5-peak 
in the range \i + -E'o(^ll) — w p i < u> < fj, + E (k\\) + u p i in favor of the growing weight of 
the hump. In the limit t± —>■ 0, i. e. u p i —>■ 0, these quasi-particles disappear and all their 
spectral weight transfers to the hump. The vanishing of the quasi-particle weight in the 
range fi + E (k\\) — uj p i < uj < fi + E (kn) + uj p i as t± — > is visible in the dependence 
of Z{k\\) on t± for k\\ = k F as shown by Eq. [29] and in Fig. [3j We thus come to the 
spectral function for t± = which has no low energy quasi-particle. In other words, the 
cross-over from the t± ^ Fermi liquid regime to the t± = non-Fermi liquid regime takes 
place through the decrease of the quasi-particle weight by closing the optical gap in the long 
wavelength plasmon mode. 

We note that numerically obtained spectral function shown in Fig. H] fulfils excellently 
the sum rule 

/oo 
A(k h uj)duj = 1, (33) 
-oo 
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(a) (b) 

FIG. 4: Spectral function A(k\\,u;) for small (u> p i = 0.2QeV) (a) and large (u p i = O.Q3eV) (b) value 
of the transverse plasmon frequency uj p i in the case kp = it /2b. Broad maxima for different values 
of the wave number k\\ follow from Eq. (|3ip . while 5-peaks are represented by their weight Z(k\\) 
according to Eq. 



with the agreement up to 10~ 4 in the whole range of the wave vector k\\, and for all considered 
values of t±. Finally, we notice that, in contrast to the quasi-particles in the range /i + 
E (kn)—ujpi < uj < fj,+E (k\\)+ujpi, the quasi-particles in the energy range \u — p, — E (k\\) \ > 
Q p i are not critically sensitive to the plasmon optical gap u p i and keep a finite intensity in 
the limit t± — > as was already shown in Ref. [6|. 

As was already mentioned in the Introduction, the main property of the above spectral 
function, namely the quasi-particles at low energies coexisting with the wide structure orig- 
inating from the collective plasmon branch, resembles to the result obtained in the early 
investigation of the isotropic ,,jellium" model within the GqWo approach by Hedin and 
Lundqvist 17, Q]- They showed that due to the finite long- wavelength minimum in the 
optical plasmon dispersion, Q p i, a quasi-particle with reduced weight appears in the region 
ji — Q p i < uj < n + Q p i, while the rest of the spectral weight is widely distributed at energies 
outside this range. 
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As was already argued in Ref. Q], the non- Fermi liquid regime for t± = is in the 
qualitative agreement with the ARPES spectra of Bechgaard salts which apparently do 
not show low energy quasi-particles 0, 0, 5|. On the other hand, the present results 
for the spectral function of the quasi-one-dimensional metal in the t± ^ Fermi liquid 
regime suggest that in (TMTSF) 2 PF 6 (for which t± = O.OUbeV and t = 0.125eV) the 
quasi-particle 5-peak with the weight of the order of 20% of the total spectral weight for 
a given value of fey is expected in the low energy range, at an energy distance of the order 
of ujpi = O.lSeV^ from the lower edge of the wide hump. A more directed experimental 
search, supported by improved energy and intensity resolutions, is very probably necessary 
for finding peaks with so weak intensities. 



Finally, we refer to the work 



26] devoted to the quasi-two-dimensional metals with the 



finite transverse transfer integral t± between metallic planes, with the main result analogous 
to ours. Namely the spectral function in this case also consists of the suppressed quasi- 
particle peak and a broad feature. Again, the RPA screened Coulomb interaction gives a 
strongly anisotropic plasmon branch dispersion of the form (115j) containing small transverse 
plasmon frequency compared with the longitudinal one. This result is in agreement with the 
ARPES spectra of quasi-two-dimensional high-T c superconductors in the normal conducting 
phase [27]. 



IV. DENSITY OF STATES AND MOMENTUM DISTRIBUTION FUNCTION 

Integrating numerically the spectral density A(ku,u) in terms of fen, we get the density 
of states for band electrons, 

1 — 

n(w) = 2k~ F So A ( k h") dk h ( 34 ) 

shown in Fig. [5] for two values of interchain transfer integral, t± = 0.025ey and 0.06eV. 
Three distinctive step singularities in n(u) originate from the edges of the corresponding 
quasi-particle 5-peak dispersions. In particular, the density of states falls from a maximum at 
the lowest energy of the fey-dependent quasi-particle 5-peak in the range u> < fi + Eo(k\\) — fl p i 
to a local minimum. Then it rises until the step discontinuity at the highest energy of the 
quasi-particle <5-peak in the energy range u < fi + Eo(k\\) — Q p i. Further on, n{uj) varies 
slowly from this discontinuity until the next one at the lowest energy of the quasi-particle 



15 




(a) (b) 
FIG. 5: Density of states n(u) for t± equal 0.025eV (a) and O.OGeV (b). 

5-peak in the energy range \i + E (k\\) — u p i < uj < fi + E (k\\) + u p i, accumulating the 
contribution from the spectral density hump in this range. Increasing further the energy 
above the third step discontinuity one comes to the minimum of n(u) at u — [i, the latter 
bearing the contribution from the quasi-particle at the chemical potential in the spectral 
function. 

The momentum distribution function 

n(k\\) = / A(k\\,uj)du> (35) 

is also calculated numerically, and shown in Figj6]for t± = 0.025eV (a) and O.OQeV (b). The 
deviation of areas below the curves (a) and (b) from the exact number of particles is smaller 
than 0.1%, indicating the highly satisfying self-consistency of the GqWq approximation. The 
momentum distribution has a qualitative behavior of the dressed Fermi liquid. It decreases 
from the maximal value at h\ = towards the step discontinuity at the Fermi wave number 
k\\ = kp. The height of this discontinuity is equal to the spectral weight Z{kp) of the quasi- 
particle 5-peak at u = fi. Fig. [6] again shows that this height decreases as tj_ decreases. 
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(a) (b) 

FIG. 6: Momentum distribution function for fcp = |r and t± equal 0.025ey (a) and 0.06eV (b) 
showing the discontinuity at hp. 

V. CONCLUSION 

The aim of the present analysis is twofold. 

Firstly, we investigate the cross-over from the specific spectral function of one-dimensional 
conducting band to that of standard isotropic three-dimensional Fermi liquid. We show that 
the absence of quasi-particle peaks is limited to the band with the strictly one-dimensional 
flat Fermi surface. Quasi-particle peaks appear immediately with introducing a finite corru- 
gation of Fermi surface, measured by finite t± in our approach. The spectral weight of these 
5-peaks for kn — kp is given by the expression (liZUl) and shown in Fig. [3J It has a non-power 
law dependence on the transverse bandwidth [Z ~ — (lntjj -1 ] in the limit t± — > 0. The 
rest of the spectral weight is carried by the wide feature in the energy range characterized 
by the plasmon energy Q p i. As it is shown in Section II, this result is to a great part ob- 
tained analytically after few technical simplifications which are well justified in the limit 
t± <C t , Vtpi. 

Although, due to this limitation, our method of calculation cannot be extended towards 
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pure three-dimensional regime (t± ~ to), the plausible expectation is that the quasi-particle 
spectral weight will increase continuously as t± further increases, approaching the three- 



dimensiona 
Lundquist 



re gim e 



13, 
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with quantitative properties obtained long time ago by Hedin and 



191 ] . It is worthwhile to stress again that, as the above Z vs t± de- 



pendence illustrates, the present calculations, unlike some others (e. g. Refs.[20i. 121]). are 
not simple power low expansion in terms of t±, and in this respect are complementary to 



the higher- dimensional bosonization approach developed in Refs. [131 . [14|. The essential 
reason for the inadequateness of the perturbation approach in terms of t±, even in the limit 
t± — > 0, is to be recognized in a qualitative change of the plasmon spectrum, namely in the 
opening of the gap in its long- wavelength limit. This gap in turn enables the appearance 
of quasi-particles in A(k, ui) already within the Go Wo approximation. The word of warning 
here concerns the applicability of the GqWq approximation itself. Strictly, it is limited to 
the range of weak screened Coulomb interaction, the relevant criterion being Q p i < t . In 
some of illustrations presented here we allow for values of Q p i above this range, expecting 
that no qualitatively new situation takes place in the intermediate range Q p i ~ to- This 
range, as well as the range of strong long range Coulomb interaction (even after the RPA 
screening taken into account) however still awaits a better understanding. 

Present analysis can also provide some estimations on the possible observability of simul- 
taneous appearance of quasi-particles and wide humps in experiments measuring spectral 
properties. The energy resolution in reported photoemission measurements on Bechgaard 



salts varied between lOmeV and 30meV 



3|, |4], |5j . Additional complication comes from in- 



dications that surface effects could have affected low energy parts of existing ARPES data 



281 ] . Thus in order to observe a dispersing sharp low-energy quasi-particle with the narrow 



width ranging up to lOmeV, it will be necessary to have an increased energy resolution 
at low energies and an enhanced bulk sensitivity of the ARPES spectra. We believe that 
such demands are achievable, particularly because our estimations suggest that the spec- 
tral weights of quasi-particle peaks are expected to range up to 20% of the total spectral 
weight, and to be positioned at binding energies ranging up to the energies of the order of 
250meV, appearing in the coexistence with characteristic wide humps already observed at 
higher energies. 

Among quasi-one-dimensional materials investigated in photoemission measurements the 
acceptor-donor chain compound TTF-TCNQ appears to be a particularly interesting ex- 
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amp 



29 



There are various indications, like e. g. the infrared optical measurements 
30, |3jJ, that it has a soft longitudinal mode at lOmeV in the metallic phase. This 



mode was explained theoretically within the model of the quasi-one- dimensional metal with 
two bands per donor and acceptor chains and the three-dimensional RPA screened electron- 
electron interaction [32|. It was shown that the appearance of such mode in the low energy 
range is due to the strong coupling between the plasmon and the collective inter-band dipolar 
mode. As for the ARPES spectra, they show the absence of the low energy quasi-particles 
and the one-dimensional dispersion of electron bands l|, [2J. However the bandwidth values 
from these data are two to four times larger than the values obtained by earlier theoretical 
and experimental estimates 33j. This signalizes that it is necessary to include electron- 
electron interactions in order to improve quantitative interpretation of the data. More 
precisely, it remains to investigate the influence of the elsewhere observed low energy mode 
on the low energy spectral properties of the quasi-one-dimensional metal with one electron 
band per donor and acceptor chains within the GqWq approximation, but with the RPA 
screened Coulomb electron-electron interaction obtained for the model with two bands per 

n n 

chain [32]. Taking into account the results we obtained in Ref. [6[ and in the present work, 
we expect that this low energy mode is also responsible for the low energy spectral properties 
of TTF-TCNQ. The full analysis of this question is under way. 
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